Domain Wall Resistance based on Landauer's Formula 
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The scattering of the electron by a domain wall in a nano-wire is calculated perturbatively to 
the lowest order. The resistance is calculated by use of Landauer's formula. The result is shown to 
agree with the result of the linear response theory if the equilibrium is assumed in the four-terminal 
£N| i case. 
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I. INTRODUCTION 



The electronic transport properties in nano-size ferromagnetic structures has been drawing growing interest recently. 
Of most interest would be the magnetoresistance effect, namely, the change of the resistance by application of a 
magnetic field. The magnetoresistance in nano-structures is due to the rearrangement of the magnetization, namely 
the creation and motion of domain walls, driven by the external field. Recently high-sensitivity measurements on 
magnetic nano-structures such as wires and contacts has become possible and intensive experimental studies of the 
magnetoresistance due to a single domain wall have begun |l]-|3| [Q-|9| . These recent experiments have evoked a renewal 
of interest in theories of the scattering of the electron by a domain wall fl|l0|-[i"5| |]l6| , |l7f , although the first theoretical 
study of this problem goes back to more than 25 years ago fll] . 

The scattering mechanism is due to the non-adiabatic deviation of the electron spin from the local magnetization 
as the electron goes through the wall. In the case of a wall with thickness A S> kp l (hp being the Fermi wave length), 
the electron spin can adiabatically follow the spatial rotation of the magnetization, resulting in a negligiblly small 
reflection. This is the reason the wall contributes very little in the magnetoresistance in the 3d metals in the bulk 

^: » 

The situation changes in nano-magnets, where the profile of the wall can be different from that in the bulk magnets. 
Q ■ The extreme case would be an atomic size contact. There the wall is believed to be constrained in a small contact 
region of an atomic scale due to a rapid change of the contact shape [p^[ , and thus the adiabaticity does not hold 
resulting in a strong scattering Jl4|,|7j . Indeed an experiment on Ni contact revealed a magnetoresistance of about 300% 
I ' at room temperature ||, which is significantly large compared with other giant magnetoresistive materials such as 
, multilayers. The result has been interpreted as due to a strong electron scattering by a domain wall, which is trapped 
in the narrow contact region in the absence of the magnetic field, and which is expelled by applying a magnetic field 

In this paper the domain wall resistance in nano-wires or contacts is calculated by use of Landauer's formula. The 
advantage of this formalism is that we can see how the density of states of the sample and the lead appears in the 
resistance, which is not possible in he linear response calculation as carried out in ref. Q. The system is treated as 
in one dimension and the calculation is at zero temperature. The transmission coefficients are calculated in section 2, 
and the resistance is calculated in the case of four- and two-terminals in section 3. The result is shown to agree with 
the result of the Mori formula of the linear response theory in the case of the equilibrium flow in a four-terminal case. 
Throughout the paper we consider only channels which contribute to the transport and so do not discuss the 
I ■ quantized behavior of the conductance. 

The phase of the transmission coefficient contains important information on the stability of the electron system. 
This has been considered elsewhere [^Cj]. 
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II. SCATTERING PROPERTIES 



We consider the exchange interaction between the electron spin and the local spin. The Hamiltonian is 

H = e k cl a c ka - g J dxS(x)(c< crc), (1) 

ka 

where ej, = h 2 k 2 /2m — tp (ep being the Fermi energy). The spin index is denoted by a = ± and er is the Pauli 
matrix (spin indices are suppressed in the second term) . The local spin S has a spatial dependence of a domain wall. 
In terms of a polar coordinate (S z = ScosB) the wall is represented as cos 6* — > 1 for x <C — A and cos — > —1 for 
x ^> A, A being the thickness of the wall. We consider two typical profiles of the wall; 
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bulk wall : cos 9 
linear wall : cos 8 



tanh(a;/A) 

1 (x > A) 
x/X (-A < x < A) 
-1 0<-A) 



(2) 
(3) 



where the first case corresponds to the wall determined by the uniaxial anisotropy energy, which would be the case 
of bulk magnets, and the second linear profile would be realized in small contacts. 

To proceed we carry out a local gauge transformation in the electron spin space so that the z-axis is chosen to be 
along the local direction of spin S pl| . The transformation is written as a a — (Uc) a , where U is a position-dependent 
matrix defined as U = cos (|) a z + sin (|) a y and the electron in the new frame is denoted by a a . The Hamiltonian 



is transformed to be H = Hq + _ff; nt where Hq 



e k — cAis the energy of uniformly polarised 



tkc,a k(J a k <j and e ka 

ka 

electron with the exchange splitting A = g |S|. The term H lnt describes the interaction between the electron and the 
wall and is given by O] 



Hi, 



h 
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kq 



A q a\ +q a x a k + — ^ A p A. 



where 



dxe~ tqx XIQ, 



(4) 



(5) 



is a domain-wall form factor. Note that V8 — ► for |a;| — + oo and so the interaction is localized near the wall. The 
spin-conserving process in the rotated frame as the electron goes through the wall corresponds to the adiabatic process, 
where the electron spin follows the rotation of the magnetization, and the spin flip corresponds to the deviation from 
the adiabaticity. 

The interaction Hi nt leads to a scattering of the electron and thus contributes to a resistivity due to the wall as has 
been discussed before based on a linear response theory fll]Jl5| ,[7[ . Here we take a different approach of calculating the 
T-matrix perturbatively with respect to H mt . The solution of the Schrodinger equation with the energy of E, which 
we denote by |*.e > (two component vector corresponding to two spin state), is given as 



\*e >=\^ >+G° E T E \^ >, 



(0) 



where |*o > is the incoming wave, G° E = H _ E _ i0 is the free Green function with energy E and T-operator is defined 
as T E = -Hint i u. — ■ The small imaginary part in G E is to choose the correct boundary condition of outgoing wave. 

int 

We consider the incoming wave with energy E and spin a. They are explicitly written as 



< a; 1*1 > 



ik+x 





ik— x 



(7) 



where k a = y^2m(E + crA) is the wave vector of the spin a component. The spin a' component of the full wave 
function for the incident wave > is then obtained from eq. (0) as 



*£°(x) = S^e*-* + [ dx'G° EtT ,(x - x')Y,e lk ' X 'Tk>*',k„*, 



(8) 



where G Ecr ,(x - x') =< xa'\G E \x'a' > 
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e ik a ,\x-x'\ and T k _< k ' a '\T E \k a a >. 



The transmission and reflection coefficients are defined from the asymptotic behavior as 



The result is 



Sv'ae " x + r a > a e tk °' x (x — > -oo) 
t rT , a e ik "' x (x -> +oo) 
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(9) 
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The diagonal (in spin) component of the T-matrix for the incoming wave vector of k a is calculated as 

Tk„a,k„a = < k a ,a\H int \k a ,a > + ^ < k a ,a\H- mt \k' >< k' , -a\G%\k' , -a >< k' , -a\H int \k a , a > 

k' 

+0(A 4 ) 
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1 - 



{q + k a ) 2 - k 2 _ a - iO 



+ 0(A 4 ). 



(11) 



Off-diagonal components are given by 



- ife— crj— (T^ertT ^TftL 

The result of the transmission and reflection coefficients are 



, /, , s . \k(x + k-a) . 
t-a,a{k a ) = ~l ~ A^ ka+k _ o 



f—u,cr\ka] 



4fc__ 



4fc_ 



A. 



r«AK) = 0(A 2 ). 



(12) 



(13) 



The scattering probabilities are defined as T a i a = \t<j'a\ 2 and R a '<r = |fVo-| 2 - These probabilities for the electron at 
Fermi level (i.e., k a = kFcr) are obtained to the second order in A as 



T„„ = 1 - 



4(1 -C 2 ) 



(C 2 |A + | 2 + |A_| 2 ) 



4(1 - aC) 2 
1 

4(1 - aC) : 



L4_| 2 



(14) 



and R aa = 0(A 4 ), where A+ = A kF+ + kp _ (= ^42fc F ) and A_ = Afc F+ _fc F _ (= Ai kF Q) describes a backward and forward 
scattering with spin flip, respectively ( hp = (&f+ + kp-)/2 and £ = (fc_F+ — kp-)/(kp + + kp-)). The conservation 
of the current, 



^Fcr — kp a T a<7 + kp.—o-T—o-fj + kp— a R— ar7 , 



(15) 



is satisfied. 



III. WALL RESISTANCE 



To relate these quantities to resistance, we consider a four-terminal case where the system is coupled with two leads, 
L\ and Li and two terminals, Ta and Tb (Fig. 0). (The gauge transformation does not affect the resistance, since 
the shift of the current due to the transformation vanishes far away from the wall.) L\ and Li are connected with a 
reserver (current source) and Ta and Tb are to measure the voltage difference across the sample. The argument goes 
similarly to the non-magnetic case described in refs. p2| , p3| but taking account of spin dependence. The chemical 
potential of L\ and Li are fixed to be \x\ = tp + eVb and fi2 = ep (Vq is the applied voltage), respectively. The 
chemical potential of Ta and Tb, ha and Hb are determined below by the condition of vanishing current flow through 
each terminal. The true voltage difference across the sample is different from Vo and is given as V = (fiA — ^b)/z- 
The resistance of the system is written as i? w = V/I, I being the current through the sample. 

We first consider the case where the leads are ferromagnets with density of states for the spin state a given by 
N a . The case of non-magnetic leads is discussed at the end of this section. The number of the electron with spin a 
injected from L\ is N a eVo, and thus the electric current is written by use of transmission probability as 
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I = eVY J Na-^j j k F ^T a , a . (16) 

The voltage ha and \xb are determined by the balance between the terminal and lead. In L\, the density (per energy) 
of spin a electron is written as 

U Xa = \{N„{1 + Raa) + N^R^-a). (17) 

Here a factor of 1/2 comes from choosing right moving electron. The density at the other lead is given by 

n 2 * = ^{N a T aa + N_„T^ a ). (18) 
The condition of vanishing flow between L\ and Ta and L 2 and Tg are given as 

^ ^ia(Mi - Va) = ^2i N <y - ni„){nA - M2) 

y^n 2a {Hi - Mb) = ~ n 2 CT )(Ms - M2), (19) 

respectively. This leads to the voltage difference across the sample of 

V =- y^(n la - 7i 2ct )(mi - H2) 

By use of eq. ( |l4|) the resistance to the lowest order in domain wall scattering is obtained as 

ml 1 1 

R - = ^ z„^(i-<¥ e{1+(m ^ - {e + mA ~ n - (21) 

where 5 = (N+ — N-)/(N+ + N-). Naively it appears correct to choose N a = ^jfi™ = D a , equal to the one- 
dimensional density of state. This is wrong if the equilibrium situation with a steady current flow is considered. In 
fact, the electron with down spin, which has a smaller velocity of , cannot leave the scatteringregion quickly after 
the scattering and thus the density becomes higher by a factor proportional to (see Fig. 0). Therefore under 

a steady flow, the effective density of states becomes proportional to D a /kp a , namely, N a = D a (i+aC) ' wnere 

a prefactor is chosen to conserve the total density, i.e., N+ + N- = D+ + £)_. This point has not been discussed in 
the multi-channel case carried out so far Q • Then S = — and we obtain 

R^ = ^^^r ) U + \ 2 + \A-\ 2 ]- (22) 

This result is the same as that obtained based on a linear response theory where the equilibrium is assumed. 

In the case of two-terminals, the system is out of equilibrium and the resistance is modified. The conductance in 
this case is given as the transmitted current, eq. ([If]) with N a = D a , divided by the applied voltage, V. Thus 

(7 ' 

The wall resistance is written as = G 1 — G^ 1 , where Go = ^- is the conductance without the wall. Thus to 

the lowest order in the wall scattering, 



This resistance contains only backward scattering kp a — » —kp^ a . This result is obtained from the four terminal case, 
eq. ©, by putting N a = D a . 
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Two resistances, ( |22[ ) and ( J24j ) (dashed and solid lines, respectively), are plotted in Fig. |3| for the two profiles, (|2j) 
and (H). In the bulk case, @), A q is given by 

7T 

Aq = cosh(7rgA/2) ' (25) 
and so resistance decays rapidly for \q\ > A" 1 . In the case of a linear wall, (||), decay is slower (Fig. ||); 

sin qX 

A =7 ^f. (26) 

This difference in the form factor leads to a slightly different resistance, but strong suppression of the resistance for 
kp\ 1 is universal. The difference between the two- and four-terminal cases is small if the Zeeman splitting is large 
but becomes manifest when £ is small. 

Let us mention the case where the leads connected to the sample are non-magnetic metals. If we write the hopping 
probability between the lead and the sample as , the balance of the density of states requires DT^ = Ng, where 
D is the density of states in the lead. Thus the voltage across the sample is determined by the same eqs. as (|l7|)-(p0|), 
and so the resistance is equal to that for the ferromagnetic lead discussed above. 



IV. SUMMARY 



The scattering of the electron by a domain wall is studied theoretically. By use of a local gauge transformation, 
the interaction between the electron and the wall is treated perturbatively to the lowest order. The transmission 
coefficients are obtained and the resistance was calculated using the Laundauer's formula with four-terminals. The 
comparison with the result of linear response theory was made. 
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FIG. 1. The system coupled with four terminals. L\ and Li are leads with chemical potential \i\ = €f + eVb and /12 = £f, 
respectively, which are coupled with an infinitely large current source. Ta and Tg are the terminals to measure voltage whose 
chemical potential is denoted by \xa and Hb, respectively. The voltage of the sample is given by V = (/ia — Ms)/e. 

FIG. 2. Schematic picture of an equilibrium flow in the scattering of particles with two components with different velocities. 
If some portion of the incoming particle with velocity hkF+/m is scattered to have a smaller velocity of hkir—/m, the density 
of these outgoing particle becomes higher. 

FIG. 3. Resistances by the wall for two different profiles of the wall (bulk wall given by eq. (0) (upper graph) and linear 
wall, (0) (lower graph)) plotted as functions of A for £ = 0.3 and 0.8. Solid and dashed lines corresponds to the four-terminal 
case (R„) and two-terminal case (i?„ ), respectively. 
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